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00 . To my teacher S.P.Novikov on his 75th birthday 

In this paper we construct examples of commuting ordinary scalar differential operators 
with polynomial coefficients that are related to a spectral curve of an arbitrary genus 
. g > and to an arbitrary rank r > 1 of the vector bundle of common eigenf unctions of 

the commuting operators over the spectral curve. This solves completely the well-known 
existence problem for commuting operators of arbitrary genus and arbitrary rank with 
polynomial coefficients. The constructed commuting operators of arbitrary rank r > 1 and 
arbitrary genus g > are given explicitly, they are generated by the Chebyshev polynomials 

m '■ T r(x). 

The first operator L 2r from the constructed commuting pair of genus g and rank r is of 
^ \ order 2r and has the form 

where T r (x) is the Chebyshev polynomial of the first kind, of degree r, r > 1, the notation 
^ \ T r (^J means the ordinary differential operator, which is the Chebyshev polynomial T r of 

£, a is an arbitrary nonzero constant, b is an arbitrary constant. The second operator 
M(2 9 +i) r from the constructed commuting pair of genus g and rank r is of order (2g + l)r 
and there exists a certain hyperelliptic relation 

M(2g+i) r = L^ r +l + a2 9 L^ r + . . . + aiL 2r + a , (2) 



<(2g+l)r 

where a,- are some constants 



[L 2r , M(2g+l)r] = 0, 



the coefficients of the operator M^g+iy are expressed polynomially in terms of the coeffi- 
cients of the operator L 2r and their derivatives. 



1 This research was supported by the Russian Foundation for Basic Research (grant nos. 11-01-00197 and 
ll-01-12067-ofi-m-2011) and by the "Leading Scientific Schools" programme (grant no. NSh-4995-2012.1). 
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To every generic point P = (A,/i) of the hyperelliptic spectral curve 

^2 = X 2g+1 + a ^ x 2g + _ _ _ + ^ + Qq 



(3) 



of genus g there corresponds r-dimensional space of common eigenfunctions ip(x,\,fi) of 
the constructed commuting operators L 2r and M^ g +i) r of genus g and rank r, 



Recall that up to now no examples of commuting ordinary scalar differential operators 
that are related to a spectral curve of genus g > 1 for rank of the form r = 6s ± 1, s > 1, 
were known. For all other values of genus g and rank r, explicit examples of commuting 
operators even with polynomial coefficients were constructed (see [1] and references therein 
and here in the text below). We conjectured in [1] that there exist commuting ordinary 
scalar differential operators with polynomial coefficients that are related to a spectral curve 
of an arbitrary genus g > 1 also for an arbitrary rank of the form r = 6s ± 1, s > 1, and in 
this paper we construct such examples using recent Andrey Mironov's remarkable results 
on self- adjoint commuting ordinary scalar differential operators of rank 2 and an arbitrary 
genus g (see [2], [3]). Thus, examples of commuting ordinary scalar differential operators 
with polynomial coefficients are constructed now for arbitrary genus g > and arbitrary 
rank r > 1. 

The study of the commutation equation for two scalar ordinary differential operators 
is one of the old classical problems of the theory of ordinary differential equations (see 
[4]-[8]). The problem of finding the general form of commuting operators of fixed genus 
g and fixed rank r, and even constructing partial examples of such commuting operators, 
is very nontrivial for g > and r > 1. Moreover, there was the very interesting problem 
of existing commuting ordinary scalar differential operators with polynomial coefficients 
for fixed genus g and fixed rank r (in particular, this existence question was posed by 
I.M.Gelfand at the famous Gelfand seminar at Moscow State University in 1981). It is well 
known that commuting ordinary scalar differential operators with polynomial coefficients 
give commutative subalgebras of the Weyl algebra W±, i.e., the algebra with two generators 
p and q and the relation \p,q] = 1, so they are of a special interest and we also give them 
a special attention. These problems are very well known and we will give here only very 
brief summary of the earlier obtained results. 

Let us consider a system of nonlinear ordinary differential equations on the coefficients 
of two scalar ordinary differential operators 



L 2r 1p = ty, M {2g+ l) r il = flip. 




n 



rn 



(4) 



that is equivalent to the commuting condition 

[L, M] = 0. 
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The operators are usually assumed to be in the standard canonical form (Burchnall, 
Chaundy, [6]), i.e. 

u n (x) = v m (x) = 1, u n --i_(x) = 0. (5) 

For a pair of commuting operators this can always be achieved by a change of variables 
and a suitable conjugacy (Burchnall, Chaundy, [6]). 

By the Burchnall-Chaundy lemma [5] , [6] any pair of commuting ordinary scalar differ- 
ential operators L and M is connected by a certain polynomial relation Q(L, M) = given 
by the spectral curve Q(\, /i) — of the pair of commuting operators: Lip = Mip = /iip, 
and common eigenf unctions of the commuting operators define an r-dimensional vector 
bundle over the spectral curve (the spectral curve T defined by the relation Q(X,fi) = 
is irreducible and is completed at infinity with a unique point P ; the dimension r of the 
vector bundle of common eigenfunctions of a pair of commuting operators over the spec- 
tral curve at a generic point of the curve T is called the rank of the commuting pair of 
operators; the rank of any pair of commuting operators is a common divisor of the orders 
of these commuting operators). For commuting operators of relatively prime orders the 
rank is equal to 1. In the case of rank 1 the commutation equation has been integrated in 
[5]-[8]. The common eigenfunctions and the coefficients of commuting operators of rank 1 
are expressed explicitly in terms of the theta-function of the spectral curve [9]. The case 
I > 1 for spectral curves of nontrivial genus g > is much more complicated and much 
more interesting. 

The first examples of commuting ordinary scalar differential operators of the nontrivial 
ranks 2 and 3 and the nontrivial genus g — 1 were constructed by Dixmier [10] for the 
nonsingular elliptic spectral curve /i 2 = A 3 — a, where a is an arbitrary nonzero constant: 





+ X s + a] + 2x, (6) 



d \ q \ i d \ „ d 



M= i \di) + a ) + 31 UJ + 3 ^ + 3x(x + a) ' (7) 

where L and M is the commuting pair of the Dixmier operators of rank 2 and genus 1, 

M 2 = L 3 - a, [L, M] = 0, 
the orders of the commuting operators L and M are 4 and 6, rank 2; 
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M =[[iz) +x 2 + a) + + 3(x 2 + a)-^- + 3x, (9) 

where L and M is the commuting pair of the Dixmier operators of rank 3 and genus 1, 

M 2 = L 3 — a, [L, M] = 0, 

the orders of the commuting operators L and M are 6 and 9, rank 3. 

These remarkable unusual examples were found by Dixmier as commutative subalgebras 
of the Weyl algebra W\ by a quite nontrivial, purely algebraic way without any connection 
to the spectral theory of commuting operators [10]. Both the Dixmier examples are in the 
standard canonical form. 

The general classification of commuting ordinary scalar differential operators of non- 
trivial ranks r > 1 was obtained by Krichever [11]: a pair of commuting operators of rank 
r is determined by specifying the curve T, a point P e T, a local parameter k^ 1 (P) in a 
neighbourhood of P , by specifying r 2 g (g is the genus of the curve T) constants (0:^,7,), 
1 < i < rg, < j < r — 2, called the Tyurin parameters, and by specifying r — 1 arbitrary 
functions Wj(x). 

The common eigenfunctions ipj(x, P; Xq), < j < r — 1, P 6 T, normalized by the 
condition 

{^x i ^ j(yX,P]X °^) =6ij ' °-^- r_1 ' 

have the following analytic properties (Krichever, [11]): 

1. The common eigenfunctions ipj{x, P; xq), < j < r — 1, are meromorphic on the 
spectral curve T outside P , and each one has rg simple poles 7i(xo), with 

tl)ij{x,X Q ) . . 

yj{x,z;x Q ) p-r, 0<j<r-l, 1 < « < rg, 

in a neighbourhood of the pole 7«(xo). 

2. All the residues ipij(x,xo) are proportional to one of them: 

ll)ij(x,X Q ) = (Xijixotyir-^XjXo), < j < T - 2. 

3. If k~ l (P) is a local parameter on V in a neighbourhood of P , then we have the 
asymptotics 

if(x,P;x ) = |^£(V)Ar s j %(x,k;x ), 

— » — * — * — # 

where tp(x, P; x ) = (^ (x, P;x ),..., ip r -i(x, P\ x Q )), £ (x) = (1, 0, . . . , 0), £ s (x ) = for 
s > 1, and &o(x, k; x ) = (3>o ) is a solution of the matrix equation 



dx 



= S$ , &j(x ,k;x ) = 5^, 0<i,j<r-l, 



and 
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r-2(x) 





/ 





\k + w (x) wi(x) w 2 (x) 

The analytic properties 1, 2 and 3, the arbitrary constants (7^,0;^), 1 < i < rg, < j < 
r — 2, and the arbitrary functions w (x), . . . ,w r _ 2 (x) determine a vector-valued function 
ip(x, P; x ) and a commuting pair L, M of rank r and genus g in general position (Krichever, 
[11])- 

In order to find commuting operators Krichever and Novikov proposed in [11], [12] the 
method of deforming the Tyurin parameters (ji(xo), aij(xo)) that allowed to obtain the 
general form of commuting ordinary scalar differential operators of rank r = 2 for arbitrary 
elliptic spectral curve (genus g = 1) [13]. 



Let us consider the Wronskian matrix 



ij(x,P;x ) 



( ^0 

^0 



Ik 



V 1>< 



(r-1) „/,(r-l) 



— # 

of the vector- valued function ip{x, P; x ), and 



/ 





^ 1 



1 






1 



••• 1 

\ XO Xl X2 ■■■ Xr-1 J 



4-^ 



\ 




/ 



where \j{ x )P) are meromorphic functions on V. 

For X — Xq the poles of Xji x o,P) coincide with 71 (x ), . . . , 7 rg (x ), and the ratios of 
the residues of the functions Xj( x o,P) at the points ji(xo) coincide with the parameters 
aij(xo): 



Ctij{x Q ) 



reS 7i(^o)Xr-l 

In a neighbourhood of Pq on V the functions Xj( x o, P) have the form 

X o(x,P) = k + w {x) + O{k- 1 ), 
Xs{x,P) =w s {x) + 0{k- 1 ), l<s<r-2, 
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Xj (x, k) = C ^ [X + dij{x) + 0(k - 7i (x)), 





< J ' < r — 1, 1 < i < rg. 



Krichever and Novikov proved [11], [14] that the functional parameters 7i(x) and ctij(x), 
1 < % < rg, < j < r — 2, satisfy the system (the Krichever-Novikov system of equations 
of deformation of the Tyurin parameters) 



For arbitrary functions Wi(x), < i < r — 2, the solution of the Krichever-Novikov system 
uniquely determine an array of meromorphic functions Xji x o, P) with the required analytic 
properties (Krichever, Novikov [11], [14]). 

By the asymptotics of Xj( x o^P) i n a neighbourhood of P on T one can reconstruct 
the commuting operators. We note that the coefficients of one from a pair of commuting 
operators are always expressed polynomially in terms of the coefficients of the second 
operator of the commuting pair and their derivatives. 

The general form of commuting ordinary scalar differential operators of rank 2 for 
an arbitrary elliptic spectral curve was found by Krichever and Novikov [13]. The func- 
tional parameter corresponding to the Dixmier example of rank 2, genus 1 among all the 
Krichever-Novikov commuting operators of rank 2, genus 1 was found by Grinevich [15]. 
The class of self-adjoint Krichever-Novikov commuting operators of rank 2, genus 1 was 
singled out by Grinevich and Novikov in [16]. We found explicit formulae for self- adjoint 
commuting operators of rank 2, genus 2 depending on two functional parameters connected 
by a third-order nonlinear ordinary differential relation (see [17]), but this relation was not 
resolved. The general form of commuting ordinary scalar differential operators of rank 3 
for an arbitrary elliptic spectral curve (the general commuting operators of rank 3, genus 
1 are parametrized by two arbitrary functions) was found by Mokhov [17], [18], where 
also the functional parameters corresponding to the Dixmier example of rank 3, genus 1 
among all the commuting operators of rank 3, genus 1 were found and all commuting op- 
erators of rank 3, genus 1 with rational coefficients were singled out. Moreover, examples 
of commuting ordinary scalar differential operators of genus 1 with polynomial coefficients 
were constructed for any rank r (some commutative subalgebras of the Weyl algebra W\) 
(see also [19] and other related papers [25]-[33]). The description of commuting ordinary 



1i Qr— I; 



a i0 ~ a iO a ir-2 + ^iO^ir-l 



dij, j > 1. 
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scalar differential operators with polynomial coefficients (commutative subalgebras of the 
Weyl algebra W\) is a separate nontrivial problem, and this problem has not been solved 
completely yet even for the Krichever-Novikov commuting operators of rank 2, genus 1, 
which are rationally parametrized by one arbitrary function (this problem was considered 
and studied in [20], [21]). 

Recently Mironov [2] (see also earlier papers [22]-[24]) constructed for any genus g > 1 
remarkable examples of commuting ordinary scalar differential operators of ranks 2 and 3 
with polynomial coefficients that generalize naturally the Dixmier examples of ranks 2 and 
3, genus 1: 

L =\(t-) +x 3 + a) + g(g + l)x, (10) 



M 2 = L 29+1 + a 2g L 29 + ... + ai L + a , (11) 

where Oj are some constants, a is an arbitrary nonzero constant, L and M are the Mironov 
commuting operators of rank 2, genus g (the orders of the operators L and M are 4 and 
4g + 2, respectively), the coefficients of the operator M are expressed polynomially in terms 
of the coefficients of the operator L and their derivatives, [L, M] = 0; 

^((£) 3+ * 2+ «) 2+ ^ +i ^' < i2 > 

M 2 = L 29+l + a 2g L 29 + . . . + a x L + a , (13) 

where a« are some constants, a is an arbitrary nonzero constant, L and M are the Mironov 
commuting operators of rank 3, genus g (the orders of the operators L and M are 6 and 
6(7 + 3, respectively), the coefficients of the operator M are expressed polynomially in terms 
of the coefficients of the operator L and their derivatives, [L, M] = 0. 

Using Mironov's results, in our paper [1] we constructed examples of commuting ordi- 
nary scalar differential operators with polynomial coefficients that are related to a spectral 
curve of an arbitrary genus g for an arbitrary even rank r = 2k, k > 1, and for an arbitrary 
rank of the form r = 3k, k > 1. 

The operators L and M of orders 4k and 4kg + 2k, respectively, 

2 

+ 9 ( 9 + l)|, (U) 

M 2 = L 29+1 + a 2g L 29 + . . . + ai L + a , (15) 

where Oj are some constants, a is an arbitrary nonzero constant, are commuting operators 
of rank r = 2k, k > 1, genus g, [L, M] = 0, the coefficients of the operator M are expressed 
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polynomially in terms of the coefficients of the operator L and their derivatives [1]. For 
k > 2 the commuting operators L and M have the standard canonical form. 
The operators L and M of orders 6k and 6kg + 3k, respectively, 

- Kir - Kir + &2 G§r 1 + 

k - 2 





+k{k ~ l) \di) +(-) -■'' :i + n I ~^ Dr. (1C.) 

M 2 = L 29+1 + a 29 L 29 + . . . + ai L + ao, (17) 

where dj are some constants, a is an arbitrary nonzero constant, are commuting operators 
of rank 3k, k > 1, genus g, [L,M] = 0, the coefficients of the operator M are expressed 
polynomially in terms of the coefficients of the operator L and their derivatives [1]. For 
k > 1 the commuting operators L and M have the standard canonical form. 

In [2] Mironov proved that for commuting operators L, M of rank 2 with hyperelliptic 
spectral curve T of arbitrary genus g, 

r : /i 2 = X 2g+1 + a 2g X 29 + . . . + aiA + a , (18) 

where a« are some constants, the operator L of order 4 is self-adjoint if and only if 

Xi(x,P) = Xi(x,cr(P)), 

where a is the involution 

<r(A,/i) = (A, -fi) 

on the hyperelliptic curve T. 

Any self-adjoint operator L of order 4 in the standard canonical form can be represented 

as 

L= l — + V(x)J +W(x). (19) 

By the Mironov theorem [2], if the operator L of order 4 from a pair of commuting 
operators L, M of rank 2 with hyperelliptic spectral curve T ( fl8l) of arbitrary genus g 
is self-adjoint, then the operator L, obviously, has the form ( Tl9|) and the corresponding 
functions Xo( x , P), Xi( x ^P) have the form 

^ = -k<!fa) + Q§T)- y to ^-W^y (20) 
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where P = (A,/i) and Q(x, A) is a polynomial in A of degree g with coefficients depending 



on x, 



Q(x,X) = (X-X 1 (x))---(X-X 9 (x)), (21) 
such that the following relation holds: 

(A - W(x))(Q(x, A)) 2 - V(x)(Q x ) 2 + \{Q XX ) 2 - l -Q x Q xxx + 

+Q [V X Q X + 2V(x)Q xx + ^Q xxxx ^j = \ 2g+1 + a 2g X 2a + ... + ai X + a . (22) 

These results allowed to prove the following theorem (Mironov, [3]): the self-adjoint 
operator 

( d 2 Y 
L = i— + aiS(x) + a J + a x c 2 g(g + l)S(x), (23) 

where the function S(x) satisfies the equation 

(S x ) 2 = c 2 (S(x)) 2 + Cl S(x) + c , (24) 

a\ and c 2 are arbitrary nonzero constants, ao, c\ and cq are arbitrary constants, form a 
pair of commuting operators L, M of rank 2 with a certain hyperelliptic spectral curve T 
fTT8|) of genus g. 

In particular, if S(x) = coshx, c 2 = 1, c\ = 0, c = — 1, then all the conditions are 
satisfied and the operator 

\dx 2 



L = ( — - + a\ cosh(x) + «o ) + ot\g{g + 1) coshx, (25) 



is the first operator of order 4 from the commuting pair L, M of genus g and rank 2 with 
a hyperelliptic spectral curve. 

Let us consider the change of variable 

x = lnP(z), 

where P(z) is an arbitrary nonconstant function. Then the operator (1251) is transformed 
to the form 

f (p(z)) 2 d 2 p(z)((p z ) 2 - p(z)p zz ) d (p(z)) 2 +i y 

L= l7^F^ + (pj 3 Tz + ai 2P(z) +ao ) + 

+«ig(9 + l) (f ^ ) +1 - (26) 



9 



In particular, if P(z) is an arbitrary nonconstant rational function, then we get an operator 
of order 4 from the commuting pair of genus g and rank 2 with rational coefficients, but 
this operator is not in the standard canonical form. 
Let us consider the special case 

P( z ) = ( z + Vz 2 - l) r , r = ±l,±2,... 

Then the operator fl26|) can be represented as 

L = ( (1 " z2) ^ ~ z Tz + aTr(z) + b ) ~ ar2 ^ 9 + l ) T '( z )> ( 27 ) 

where r is an arbitrary nonvanishing integer, T r (z) is the Chebyshev polynomial of the first 
kind, of degree r or — r, a is an arbitrary nonzero constant, b is an arbitrary constant. The 
operator (I2T1) is the first operator of order 4 from the commuting pair L, M of genus g 
and rank 2 with a hyperelliptic spectral curve for arbitrary nonvanishing integer r, but it 
is not in the standard canonical form. We note that all these operators have polynomial 
coefficients and give commutative subalgebras of the Weyl algebra W\. 
Recall that 

T (z) = l, T 1 {z) = z, T r {z) = 2zT r . 1 {z)-T r . 2 {z), T_ r {z) = T r {z). (28) 

It is very interesting that the Chebyshev polynomials of the first kind T r (z) are commuting 
polynomials 

T n (T m (z)) = T nm (z) = T m (T n (z)). (29) 

We consider that it is no mere chance. Earlier we noted an important role the Chebyshev 
polynomials of the first kind in the theory of commuting operators with polynomial coef- 
ficients and with elliptic spectral curve (see [21]). It would be very interesting to clarify 
these interconnections. 

After a natural automorphism of the Weyl algebra W\ we obtain the operator ([1]) from a 
commuting pair of rank r and genus g with polynomial coefficients in the standard canonical 
form. 

Theorem. The operators L 2r and M^ g +i)r °f orders 2r and (2g + l)r, respectively, 

L » = { aT - (I) " ~ Sx i + 12 + ") 2 - ar2 ^ + 1)T " (s) • (30) 

M (2g+l)r = L 2r +1 + Q 2ff^2r + • • • + a l L 2r + Ofl, (31) 

where a{ are some constants, T r (x) is the Chebyshev polynomial of the first kind, of degree r, 
r > 1, the notation T r (^) means the ordinary differential operator, which is the Chebyshev 
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polynomial T r of a is an arbitrary nonzero constant, b is an arbitrary constant, are 
commuting operators of rank r, genus g, [L 2r , M( 2 g+i)r] — 0, the coefficients of the operator 
M(2g+\) r are expressed polynomially in terms of the coefficients of the operator L 2r and 
their derivatives. For r > 3 the commuting operators L 2r and M( 2g+ i) r have the standard 
canonical form (for a = l/2 r ~ 1 ). For r = 1 this pair of operators is commuting one of rank 
2 and genus g. 

Examples. 

1) Rank 4, genus g: 



4 / , \ 2 x 2 



L *= ill) -<* 2+i >(£) - 3 4 + * 2+q 

d \ 4 / d x 



2) Rank 5, genus g: 



+ - 5 I' (32) 



(fey iVcfey x V^ry ~ \ x wjdx 



£io = I ( — ) - t ( — 1 - ^ ( — ] - ( 3x - — ) — + x 2 + a 



3) Rank 6, genus g: 



4 / ^ \ / 7 \ 2 N 2 



dx J 2 \dx J \ 16 / \dx J dx 
d \ 6 3 / d \ 4 9 / d x 



+ -5 s Mi' (34) 



4) Rank 7, genus g: 

d\ 7 7fd\ 5 7 ( d\ 3 2 fd\ 2 / 7\ d 2 

1 -i\Tx) + sUJ -* UJ -r + 64j^ +x +a 

^ +i )iis) T -Ks) B+ Ks) a -ssi- (35) 
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Abstract 



In this paper we construct examples of commuting ordinary scalar differential operators 
with polynomial coefficients that are related to a spectral curve of an arbitrary genus 
g > and to an arbitrary rank r > 1 of the vector bundle of common eigenfunctions of 
the commuting operators over the spectral curve. This solves completely the well-known 
existence problem for commuting operators of arbitrary genus and arbitrary rank with 
polynomial coefficients. The constructed commuting operators of arbitrary rank r > 1 and 
arbitrary genus g > are given explicitly, they are generated by the Chebyshev polynomials 
T r (x). 



14 



